Orbital-limiting and modulation induced by missing parity in non-centrosymmetric 

superconductors 
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We examine the depairing field Hc2 (T) of non-centrosymmetric superconductors with a spin-orbit 
couphng larger than the Zeeman energy at Hc2{0) by taking account of the mixing of spin-singlet 
and triplet pairing states due to the missing parity. When the singlet and triplet pairing components 
are mixed with an equal weight in a cubic non-centrosymmetric system, the paramagnetic depairing 
effect is significantly suppressed so that Hc2(T) approaches its orbital-limited value. A similar event 
also occurs in a quasi two-dimensional Rashba non-centrosymmetric system. The present results 
are relevant to the H-T phase diagrams of CePtaSi and the families of Li2Pd3_a;Pta;B. 

PACS numbers: 



I. INTRODUCTION 

A Cooper-pair condensate in a superconductor is 
destabilized by two kinds of field-induced mechanisms, 
the paramagnetic depairing and the orbital depairing. 
The orbital- hmited case, i.e., a situation with a depair- 
ing field Hc2{T) determined only by the increase of the 
number of vortices, conventionally occurs for spin-triplet 
superconductors with an equal-spin pairing. Recently, it 
has been noticed that the orbital-limited case occurs in 
Rashba non-centrosymmetric superconductors [l| in an 
applied field (H || c) perpendicular to the basal plane on 
the space inversion asymmetry In contrast, the para- 
magnetic depairing is so effective in Rashba superconduc- 
tors in the parallel field H _L c that, as seen in CeRhSis |3i] 
and CelrSia [1] , Hc2 (T) is strongly suppressed, and that 
the vortex state shows peculiar modulated structures Q . 
However, it should be noted that these results in the 
parallel field were obtained in the case with a purely sin- 
glet (or purely triplet) pairing. Originally, a hallmark of 
non-centrosymmetric superconductors is a missing parity 
and the resulting coexistence of spin-singlet and triplet 
pairing symmetries Jj. However, this mixing of pairing 
symmetries and a coupling between them have not been 
taken into account so far in studying Hc2{T) and the 
vortex states occurring below it. In relation to this, it 
should be noted that the strong anisotropy of Hc2 due to 
the anisotropic paramagnetic effect does not seem to be 
a common feature between Rashba superconductors: A 
strongly suppressed Hc2 in the parallel field case was re- 
alized in a couple of materials with tetragonal struc- 
ture, while the 7Jc2-curves in CePtaSi are nearly isotropic 

In this work, we study the depairing field Hc2{T) 
of non-centrosymmetric superconductors with coexisting 
singlet and triplet pairing components. Both of systems 
with the spin-orbit coupling of Rashba type in H _L c 
and those with that of cubic type will be considered here 
because a comparison between those two cases is found 
to be useful. Following previous works [2, 5, 7], a mod- 
erately large value of the spin-orbit coupling C, will be 
assumed, Max(r, ^ibH) <^ ( <^ Ep, where ^ibH is the 



Zeeman energy of a conduction electron, and Ep is the 
averaged Fermi energy. In such materials in a single pair- 
ing channel, a one-dimensional modulation [3, of the 
superconducting order parameter A and the resulting in- 
crease of Hc2 occur as a result of the small but nonva- 
nishing Q/Ep. We show below that a coupling induced 
by spatial variations of A between the coexisting singlet 
and triplet pairing channels is a much stronger origin of 
elevating iJc2-values. When the two pairing channels are 
equally important, one of the two FSs splitted due to the 
missing parity is favored for the gap formation, and this 
imbalance between the FSs leads to a more perfect disap- 
pearance of the paramagnetic depairing than that due to 
a finite C/^f In fact, under proper conditions, the 

present mechanism based on the missing parity results 
even in the orbital-limited situation where the paramag- 
netic depairing is quenched. It is argued that the present 
results should be closely related to the fact that the Hc2- 
anisotropy in Rashba superconductors significantly de- 
pends on the materials [1, 0, 0] and are also relevant 
to the ^-dependence of pairing states of Li2(Pd3_2,Pt:r)B 

Through this paper, the 77c2-enhancement due to the 
mixing of the singlet and triplet pairing channels is dis- 
cussed, for clarity, by focusing primarily on the cubic 
case. In the cubic case, the i7c2-enhancement is accom- 
panied, as in the centrosymmetric case with the Fulde- 
Ferrell (FF) state 0] , by a helical modulation parallel to 
H of the phase of A. It will be pointed out that this 
cubic case corresponds to an ideal situation of the famil- 
iar FFLO mechanism of an 7Jc2-eiihancement in which 
the paramagnetic depairing is cancelled by a modulation 
of A. In addition, it is pointed out that, in this cubic 
case, the Larkin-Ovchinnikov (LO) state with a periodic 
amplitude modulation parallel to H cannot obtain a 
gain in energy necessary for its realization. 

This paper is organized as follows. In Sec. II, the start- 
ing model and the formulation are explained, and pos- 
sible i7c2 (T)-lines in the cubic noncentrosymmetric case 
are discussed. As an observable quantity measuring the 
broken inversion symmetry in the cubic case, a transverse 
component of the local magnetization in the vortex lat- 
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tice is discussed in sec. III. In sec. IV, the iJc2 (T')-curves in 
the Rashba superconductors under a field parallel to the 
basal plane are considered for a comparison with the re- 
sults in sec. II. A summary and some comments are given 
in sec.V. 



II. MODEL AND CUBIC CASE 

We start from the following electronic Hamiltonian 
T-Cei = ^ c^^(ek(55i^s2 + (Cgk + MbH) • CT^i^s^ ) Ck,. 



k,Si,S2 

1 

V 



+p/2,a '---ki+p/2,/3 



p,ki,k2 
X C_k2+p/2,(5 Ck2+p/2,7j 



(1) 



where fisH is the Zeeman energy, Sk is the bare band 
energy, aa,p are the Pauli matrices, V is the volume, and 
gk parametrizes the spin-orbit coupling. The gauge field 
will be incorporated later at the quasi-classical level. The 
pairing interaction is represented by 



a. (3,^ S 



(ki,k2) = -i u;y(f/(ki))a/3 (Tj(k2))5 7, 



i,j=s,t 



(2) 

where Ts{k) = lay, Tf(k) = i{<7y<7fj,) ■ (gk)p, the 2 x 2 ma- 
trix Wij is positive definite, and the index s (t) implies the 
spin-singlet (triplet) component. It has been assumed in 
eq.(2) that |^| is so large that the spin component, i.e., 
d- vector, of Tf(k) is protected by the spin-orbit coupling 
[J| . By diagonalizing the quadratic term of Tiei through 
the unitary transformation Ck./j = X]a=i 2 ^/3a(^)'^k,o, 
the single particle energy close to the Fermi surface (FS) 
a (= 1, 2) is given by Ek — (— )''C|gk|, and the interaction 
Hamiltonian (the second line of eq.(l)) consistently takes 
the form 



Hi, 



V 



(3) 



P i,j = s,t 



In the cubic case, gk simply becomes k = h/kp = 
zcos^k + sin^k (2;cos(/)k -I- ysin^k) up to the lowest order 
in k, and, following other works this model will be 
used here. The transformation matrix C/(k) takes the 
form cos(0k/2) + i sin(0k/2)(sin0kO'2: — cos^kfi/)- Then, 
^fp is expressed by 



P Y "-k+p/2,a"k+p/2,a 



V 



d-\i+p/2,adk+p/2M- (4) 



In eq.(3), 0(£h) corrections expressing an interband 
pairing were neglected in writing Hint, where Eh — 
Ma,x{nBH, T)/\C\. In the Ginzburg-Landau (GL) free 



energy given below, they would lead to a correc- 
tion term of 0{e^) which is safely negligible. Then, by 
decoupling Hint in the manner — X^i j '"'u [^^*^]^^^"'^ ~^ 
E..j[(u'~^)»jA*Aj-] - (A:*^'*) + A*^^(*) -I- h.c), the re- 
sulting i^^^^ in the cubic case can be represented as a 
functional of the order parameters A^ (a = 1, 2) defined 
on the resulting two FSs, where 



-iA, + (~l)°At 
V2 



In obtaining F^'^\ one needs to use the expression 

^"^^'''^ " ie-ek + (-l)1Cgk + MBHr 
Then, using the relation 

^ 2TTNa 

2\e\ + isgne(v • H + 2(-1)''+VbH • gk) ' 



(5) 



(6) 



(7) 



where Na is the density of states on the FS a, the 
quadratic term of the GL free energy F*^"^^ is given by 

Ft'^ = / d'rl E [iiw-')ss + {w-')tt)\Aa\^-2A:KaiU)Aa] 

^ L ^ — in 



- a=l,2 



+ [( (w-^)ss - (w-%t - 2i {w-^)st )AJ A2 + C.C.; 



(8) 



where 



/■OO 

i^„(n) = \ dp/(p;T)(cos(pv • H,)), (9) 

J Or- 



2t:T 



sinh(27rrp) ' 



(10) 



Ua = U {-YQz in a field H || z, H = -iV + 2eA, 
Q — l\v\, V is the Fermi velocity vector, ( ) denotes 

the (angle) average over each FS, and the identity D^^ = 
dpexp{—pD) was used. In eq.(9), a lower cutoff pc of 
the p- integral, which is of the order of the inverse of a high 
energy cutoff lUc, was introduced. This will be needed 
even in some of the ensuing expressions. Note that, in the 
present case with the cubic spin-orbit coupling breaking 
the inversion symmetry, the paramagnetic effect appears 
only through the Q-dependence, which simply shifts the 
gauge field parallel to H in a way dependent on FS (see 
the sign factor (—1)° in Ila). 

To obtain the depairing field Hc2{T), we express the 
order parameter as 



Aa = 



:Z4q)e'«-V'o(a:,2/) 



(11) 



[9| and diagonalize ^2'^'' with respect to Z^, where 
ipo{x,y) is an Abrikosov lattice solution in the lowest 
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{n = 0) Landau level (LL). This restriction to the lowest 
LL is justified as follows: It is well known that the para- 
magnetic depairing tends to reverse roles of the lowest 
LL and higher LLs [l^, In the present cubic case, 
however, the paramagnetic effect merely plays the role 
of modulation parallel to H and is ineffective for spatial 
variations of A perpendicular to H. In this sense, this 



situation is similar to the familiar orbital-limited case, 
and our neglect of higher (n > 1) LLs in this section 
is safely valid at least as far as focusing on properties 

(c) 

occuring near Hc2- Then, after diagonalizing F2 , the 
eigenvalue determining the Hc2{T)-lme is given by 



2{Ni+N2) 



dpf{p;T){Io{p){cos{pviiq)cos{pvi\Q) + (5iVsin(pu||g)sin(pwjiQ))) 



dpf(p; T){Io(p)[SNcos(pv\\q)cos{pv\\Q) + sm{pv\\q)sin{pv\\Q)]) 



2-1 1/2 



(12) 



r 



where Further, by reexpressing eq.(12) in terms of the zero 

21 field transition temperature Tc which is determined from 

In{p) - exp(- ^ ) Ln{p'^\v^\^/{2r%)), (13) eq.(12) in i7 = case, we find that the H^2{T)-cmve is 



H 



th = (2ei/)~^/^, vj^ {v\\) is the component of v perpen- 
dicular (parallel) to H, Ln{x) is the Laguerre polynomial, 
and 



(iv, -f N2)5 . ( + K!^L!k^(!^^^ ^ 



(14) 
I 



given by 



In 



dpf{p\T) 



1 - {Io{p){cos{pv\\q)cos{pv\\Q) + (5iVsin(pt;||g)sin(/9t;||(3))) 

2n 1/2 



X {Ii){p)[5Ncos{pv\\q)cofi{pv\\Q) + sm{pv\\q)sm{pv\\Q)]) 

I 



dpfip;T) 



dp5Nf{p-T,) 



1/2 



(15) 



if q is chosen so that the highest i/- value results in, where where 
5N = {N2 - Nx)/{Nx+ N2). The coefficient Za for the 
corresponding eigenstate is given by 



Za{q) = ^(l + (-l)'^sgn(JiV)i?-i/2)/2, (16) 



R=1 + 5^(^J^ dpf{p;T){Ioip)[SNcos{pvi\q)cos{pv\iQ) + sin{pv0)sm{pv^\Q)])^ . (17) 



r 



The off-diagonal term \{w-'^)st\/{Ni + N2) is a con- the order \SN\ln{uJc/Tc). If both {w-'^)st/iNi + N2) and 
sequence of the lack of inversion symmetry and will be \(w~^)tt — {w~^)ss\/[2{Ni + N2)] are negligibly small, 
nonzero in general [l3|- We expect it to be at most of 
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eq.(15) reduces to its (5 = case 

ln(^|r) + dpf{p-T)[l-{Up) 
X co&{pv«{q + Q)))] = 



(18) 



irrespective of |i5iV|, where it was assumed that 5N < 0. 
Then, Hc2{T) is given by eq.(18) with q + Q = 0, which 
is the purely orbital-hmited one independent of the para- 
magnetic depairing. In this specific case, Ai grows on 
cooling under the condition SN < (A^i > N2), while 
A2 identically vanishes. Further, the period of the phase 
modulation is precisely 2tt/Q — 7r|v|/(/iB^f)- Since this 
Q is nothing but the relative shift of the two FSs, this 
S — limit in the cubic case can be regarded as an 
ideal case of the familiar FFLO mechanism for an Hc2- 
enhancement: In the present case, the effect of the para- 
magnetic depairing is perfectly cancelled by the phase 
modulation (i.e., a nonzero to reach the orbital- 
limited case. As far as \{w^^)st/{Ni + N2)\ will be 
nonzero, however, the growth of A2 is not negligible even 

if Wss = Wtt- 

In place of the helical variation Aa ^ e''^, we 
have also examined the alternative z-dependence, Aa ~ 
exp{i{—)°-Qz), which is motivated by the fact that the 
gradient operator acting A^ is Ua- In this case, the cor- 
responding expression to eq.(15) is given by 



ln('^ 

Tc 



dpf{p;T)[l-{l + \SN\){Ioip))] 



dpSNf{p;T,) 



1/2 



(19) 



As this expression independent of Q shows, a complete 
orbital-limiting is always realized in this case. However, 
it can be seen that, as far as (5 7^ 0, the resulting Hc2{T) 
always lies below that following from eq.(15), implying 
that, as far as Wss 7^ Wtti the orbital-limiting with no 
paramagnetic depairing is not realized. 

Now, let us discuss about the i7c2 (2^)-curves following 
from eq.(15). In Fig.l, the 5 dependence of Hc2{T) and 
of the corresponding q{T)/Q just at Hc2{T) are shown 
by setting p^h'^^^\q) / {2itT^) = 0.4 and 8N = -0.1, 

where H^\t) is the orbital- limiting field in 2D case. 
In the purely singlet or triplet case where i5 = cxd, the 
same magnitude of the energy gap is formed on the 
two FSs, and a helical phase modulation parallel to the 
field at higher temperatures is merely a consequence of 
a nonvanishing 5N [2], while the sudden appearance of 
nonzero q near t = 0.38 indicates a second order transi- 
tion into the ordinary FF state [lB|. As the lower two 
curves in Fig. 1(a) show, however, a realistic |(SA^|-value 
(~ 0.1) does not lead to a remarkable increase of Hc2- 
In contrast, effects of a (5-reduction on Hc2 and q are 
more dramatic: Even for i5 of order unity, the Hc2{T)- 
enhancement due to the singlet-triplet mixing is much 
more remarkable than that due to a finite SN , and the 




FIG. 1: (Color online) (a) Dependences of Hc2{T) curves in 
the cubic case on S and SN, where h = H/H^^^\o) and 
t = T/Tc- The two solid curves were obtained for SN = 0, 
while the remaining ones are for SN = —0.1. The parameters 
MB^/orif^(0)/(2vrTc) = 0.4 and = (207rrc)-^ were used, (b) 
The corresponding t dependence of q/Q just at Hc2{T). 



slope of Hc2 [T) shows a subtle but visible increase below 
an intermediate temperature upon cooling. As Fig. 1(b) 
shows, this increase of the //c2-slope originates from the 
corresponding increase of \q\, i.e., the phase modulation 
of A, upon cooling, and the onset temperature of the 
Igj-growth increases with decreasing 5. However, this 
enhancement of the modulation never means that the 
paramagnetic depairing becomes more effective due to 
the missing parity. To explain this, we have also exam- 
ined the corresponding Hc2{T) resulting from the LO- 

type state A^ = L'^'^ Za{ci)[e"i'Mx,y) + e-"i-^o{x,y)] 
[lO|, where (po{x,y) may be different from ipo{x,y) but 
is under the same gauge. In this case, the corresponding 
expression to eq.(15) for determining Hc2{T) is given sim- 
ply by neglecting the two sin(pti||(7)sin(pw||(5) terms there. 
Then, by noting that the LO state appears upon cooling 
at the temperature where the 0{q'^) term of Fj'^'' changes 
its sign, it is easily seen that the onset of the LO state is 
independent of 5. In contrast, in eq.(15), an additional q 
dependence appearing as a consequence of a finite 5 (see 
the r.h.s. of eq.(15)) favors the helical phase modula- 
ton. Thus, the |g|/(5-growth enhanced by decreasing 5 in 
Fig.l is peculiar to the helical (phase-modulated) state 
and does not imply an enhancement of the paramagnetic 
depairing. Rather, the singlet-triplet mixing makes the 
LO state significanty unfavorable. Further, the curves in 
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Fig.l also show that an increase of the additional helicity 
due to a nonzero 5N |2l] becomes more remarkable as the 
singlet-triplet mixing is increased. 

Through the discussion on Hc2, we have implicitly as- 
sumed the mean field superconducting transition at Hc2 
to be of second order. To check its validity, let us con- 

(c) 

sider here the corresponding quartic GL term . The 

(c) 

expression of is briefly explained as follows: Recall 
that, in the familiar centrosymmetric case with no para- 

(c) 

magnetic depairing, F^ within the lowest LL takes the 
form J dxdy\ipo{x,y)\'^{J IldpjF{pj;p)) and is positive. 
Here, p denotes a unit vector parallel to a momentum on 

(c) 

each FS. In the weak coupling approximation, F^ in the 
present cubic case is simply a sum of contributions from 
the two FSs and is given by 



F, 



dxdy\ipo{x,y)\' 



E 

a=l,2 



Hdpj 



X cosiY,PM'i-i-^rQ))FiPf,P)), (20) 

j 

which is clearly positive since q < when 6N < 0. Actu- 
ally, we have verified that F^'^'' is always positive within 
the calculations performed by us. This fact implies that 
there is no occasion that the second order transition at 
Hc2 assumed above is preempted by a discontinuous tran- 
sition ^T3| in the present cubic noncentrosymmetric case. 



III. TRANSVERSE MAGNETIZATION IN 
CUBIC CASE 

In addition to the _ffc2-liiie, an observable measure of 
the phase modulation of the vortex state in the cubic 
noncentrosymmetric case will be needed to verify its pres- 
ence in real systems. Recently, the presence of a nonva- 
nishing component mx perpendicular to H of the local 
magnetization m in the ordinary FFLO vortex lattice has 
been stressed within the gradient expansion approach [l6| 
valid for sufficiently large Maki parameters. In the phase- 
modulated FF state, this transverse magnetization mx 
occurs from a nonvanishin g pe riodic component j|| (r) of 
the current parallel to H [8l.Tl6j . In this section, our result 
of a quantity corresponding to m_L will be shown. This 
result will clarify the implication of the 7fc2- enhancement 
due to the phase modulation in sec. II further. 

The local supercurrent j is given by 



ST. 



■ / dpf{p;T) (\p\^Lpl{x,y) 
exp(ip(vi -U + v^^iq- {-irQ)))ipo{x, v))NaZt. (21) 



5A 



(5A=0 



Below, each component of j will be expressed as 
= ^ji(K)FKe'K-, 




FIG. 2: (Color online) Calculated curves of J = 
Ko,y) Ko.y)\ at Ha2{T) obtained in terms of the 
data in Fig.l. The same symbols and the type of the lines as 
in Fig.l are used here. 



where i^K is the Fourier transform oi\ipo{x,y)\'^ j'K.is the 
reciprocal lattice vector of the vortex lattice, and 

ji(K) - i8e / dppf{p;T) ^ NaZl{v^co8{pv^^iq 

■'0 a=l.2 

- {-lTQ))h{p) exp(-p(v X K) j2)cos(pv • K/2)), 
j||(K) = 8e / dppf{p;T) ^ NaZl{{-v\^)Io{p) 

X sin(p^;,[ {q - (-1)"Q)) exp(-p(v x K),/2) 

X cos(pv • K/2)). (23) 

We note that, in centrosymmetric superconductors oc- 
curing on a single Fermi surface where q — 0, and 
Na = N{0), j_L in eq.(23) reduces to the corresponding 
expression j'"' of the familiar Abrikosov lattice in the 
n = LL in the orbital-limited case with Q = 0. Then, 
the local magnetization vector is given by (h — H)/(47r) 
and is related to the current through the Maxwell equa- 
tion V X h = 47rj. Below, we will focus on estimating 
not mj_ itself but rather the normalized j|| which will 

be defined here as J = j\^{0, Ko^y)/\jji^\o, Ko_y)\, where 



Ko,y = "'^^^/(S^^^T'if ) is the magnitude of the smallest 
reciprocal lattice vectors. The magnitude and sign of 
mj^ is determined by J. 

In Fig. 2, our calculation results of the normalized 
quantity J at i7c2 (r)-curves are shown. They are ob- 
tained by combining the data in Fig.l into the above 
expressions. The curves in Fig. 2 should be fingerprints 
of features peculiar to the vortex lattice occuring at least 
near 7?j,2 in the cubic noncentrosymmetric case. As is 
(22) found by comparing Fig. 2 with Fig.l, a growth of |g| in- 
duced by the mixing of the singlet and triplet pairings 



6 



results in a reduction of |m_L|, i.e., of the paramagnetic 
depairing, while this reduction of the paramagnetic effect 
is safely negligible in (5 = oo case with no mixing between 
the pairing channels, and |m^| monotonously increases 
upon cooling reflecting an enhancement of the param- 
agnetic effect upon cooling. Thus, the nonmonotonous 
i-dependence of J seen in (5 = 0.2 and 0.6 cases at 
lower temperatures is a consequence of the competition 
between an enhancement of the paramagnetic depairing 
upon cooling and its effective reduction due to a growth 
of 



IV. DEPAIRING FIELD OF QUASI 2D RASHBA 
SUPERCONDUCTORS 



Next, let us briefly explain the corresponding results 
for a Rashba superconductor with the basal plane per- 
pendicular to £ in H II ij (i.e., a parallel field configu- 
ration). In this case, gk = (k x z)/kp^ and the unitary 



matrix i7(k) is replaced by (l-|-i(sin(/)kcry— cos0kO'a;))/\/2. 



Then, in eq. (4) needs to be replaced by i 

consequently, the off-diagonal element Wst of the inter- 
action matrix appears only in the "intraband" terms 
of the GL free energy (see below). Further, follow- 
ing the purely singlet (or purely triplet) case we 
use a cylindrical FS extending and corrugating along 

z. Then, the factor |gk| in ^'p'' of eq.(4) may be re- 

placed by unity. The quadratic term F2 of the result- 
ing GL free energy, corresponding to eq.(8) in the cu- 
bic case, reduces to eq.(2) in Ref.Q in Wu, Wst 

limit. By noting that, in -^2^^ with A = Hzx, the 
gauge-invariant operator — iV + [f^z + {—)"-Q)x acts on 
Ao = (As + (-)"+! Af)/%/2 (a = 1, 2), it is convenient 
to express A^ in terms of LLs dependent on the two FSs 
in the manner A^ = Z]„>o Vn{z + {-TQr]j: x). 
Then, we have 



is) 



J 



(R) 



2V 



(24) 



where 



min(m,n) 



Wn,m{Q) = dzdx(p'^{z,x)(pn{z + 2Qrjj, x) ^ exp{-Q'^r'jj) ^ 



(-1)'"-Vn!m! 
{n-l)\{m-l)lll 



{V2QrHr+"'^^' , (25) 



and the zero field Tc is determined by 



2{N,+N2) 



dpf{p;Tc) 



■N2) 

r 



+ 



dp8Nf{p-T,) 



(w 



(N1+N2) 



1/2 



.(26) 



In this Rashba CclSG, Oil direct numerical evaluation is 
needed to examine Hc2- Nevertheless, the main result 
is already found in eq.(25) in the specific Wss — Wtt 
case where the "interband" term is absent: In eq.(25), 
the paramagnetic depairing appearing only through the 



Q-dependences in the interband terms is completely 
quenched irrespective of the w^t-value in the Wss Wu 
limit where the pairing occurs just on the FS 1 with a 
larger density of states. Then, the vortex state close to 
the resulting Hc2 is described primarily by the Abrikosov 
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FIG. 3: (Color online) (a) Examples of calculated Hc2{T) curves in the Rashba case with ^bHI^"\Q)/{2-kT,) = 0.4 in H ± c. 
The used parameters (|(w"^)tt - {w~^)ss\/{Ni + N2), SN) are (0, 0) (black solid curve), (0.4, 0) (red solid curve), (0.4, -0.1) 
(open circles), and (00, —0.1) (dashed curve). For simplicity, we have put Wst ~ 0. 



triangular lattice solution Ai with no additional modula- 
tion of the order parameter. Even in this Rashba case, a 
discontinuous normal to superconducting transition |12l | 
is expected not to occur, because, as seen above, the 
paramagnetic effect is lost in the limit where the pairing 
interactions in the singlet and triplet channels occur with 
a comparable weight, while the absence of the tempera- 
ture region with such a discontinuous transition was ver- 
ified in the opposite limit with just a single pairing state 
0. Therefore, an exotic sequence of structural transi- 
tions between different vortex lattices Q may not appear 
in systems with a small \{w~^)ss — iw~^)tt\ close to the 
orbital- limiting. Examples of Hc2 (T) curves obtained by 
directly diagonalizing eq.(25) in terms of the lower six 
(0 < n < 5) LLs are shown in Fig. 3. It seems that 
the Hc2{T) curves at least in t < 0.6 are not dependent 
much on whether still higher LLs are included or not. It 
is found that the Wgs — wtt dependence of Hc2 (T) curves 
is qualitatively similar to the 6 dependence in the cubic 
case, and that a slight but visible slope change of Hc2{T) 
near T ~ O.GSTc occurs for a smaller \5N\. Further, in 
the case with 5 — 00, i.e., with just a single pairing chan- 
nel, a clear slope change of Hc2{T) is detectable near 
t = 0.5 and is closely related to a structural transition 
between vortex lattices Q . This result will be relevant to 
a similar behavior seen in Hc2 data of CeRhSia in H _L c 
. Details of the corresponding vortex lattice structures 
will be reported elsewhere. 



V. SUMMARY AND DISCUSSION 

In the preceding sections, we have shown that, irre- 
spective of the form of the broken inversion symmetry, 
the mixing and the field-induced coupling between coex- 
isting singlet and triplet pairing states significantly sup- 
press the paramagnetic depairing effect and lead to an 
enhancement of Hc2- In the noncentrosymmetric sys- 
tems of the cubic or Rashba type, the paramagnetic ef- 
fect enters as an additional gauge field in the gradient 



n — — iV -I- 2eA acting on the order parameter A in the 
manner dependent on the FSs; 11 -I- Q for one FS and 
n — Q for the other. If either of the two splitted FSs 
is irrelevant to superconductivity, Q is trivially gauged 
away, and the paramagnetic term plays no roles of a pair 
breaking. In contrast, when both of the two FSs con- 
tribute to superconductivity, the gauge fields ±Q frus- 
trate with each other and are not cancelled by a gauge 
transformation so that the paramagnetic depairing ef- 
fectively works. In these noncentrosymmetric systems, 
either of the two FSs may become irrelevant to super- 
conductivity when both of a singlet and a triplet pairing 
channels have attractive interactions of the same order in 
magnitude, and then, the orbital-limited Hc2 is realized. 

It will be valuable here to explain relations of the 
present work with other previous ones addressing non- 
centrosymmetric superconductors in nonzero fields. In 
Ref. [H , the coexistence of a singlet and a triplet pairing 
channels was taken into account in the Pauli limit where 
i n (see the preceding paragraph) is replaced by the gra- 
dient V so that the vortices are ignored. Further, any 
modulation of A, i.e., any contribution of the gradient 
VA, was ignored there, and consequently, the strength 
of coupling between the two pairing states was measured 
only by \5N\ ~ C^jEp as in zero field case. However, 
this result is invalidated once a modulation of A is taken 
into account. Some results of a treatment in the Pauli 
limit taking account of contributions of VA were com- 
mented on in Ref. 17[ by focusing on the Rashba case. 
It seems that a divergence of iJc2 for a cylindrical FS, 
noted there [iTj . at an intermediate temperature in the 
case with both of a singlet and a triplet pairing channels 
corresponds to an orbital-limited situation found here in 
sec. IV by taking account of the vortices. However, any 
physical implication of the -ffc2-divergence and details of 
calculations leading to such results were not explained 
there [13] • The crucial point is that the coupling, induced 
by the magnetic field and a A-modulation, is present be- 
tween the two pairing channels even in — > limit. 

Finally, we discuss about relevance of the present re- 
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suits to real systems. As noted in sec. I, two Rashba su- 
perconductors, CeRhSia [3] and CelrSis \4^, show a strong 
paramagnetic effect in a parallel field, and their in-plane 
Hc2{0) values (in H _L c) are significantly suppressed 
compared with that in H || c. In contrast, the Hc2 lines 
in CePtsSi and LalrSia are nearly isotropic and show no 
sign of the paramagnetic depairing even in the parallel 
fields. In particular, CePtaSi has a large effective mass of 
the normal quasiparticles, and, in fact, the Pauli-limiting 
field Hp{0) was estimated to be much lower than Hc2{0) 
in all configurations [6]. It will be reasonable to, accord- 
ing to the present results, attribute the apparent absence 
of paramagnetic depairing in this material under a paral- 
lel field to a mixing of two pairing channels with a com- 
parable weight. It is quantitatively insufficient to regard 
the nearly isotropic Hc2 [6] a-s a consequence of a finite 
SN 0. The present view on the pairing state of CePtaSi, 
following from studies of the H-T phase diagram, is con- 
sistent with a recent proposal [T^ based on microscopic 
properties. However, it is unclear at present whether the 
ne glec t in the present work of an antiferromagnetic order 
[^■Tisj existing in CePtsSi is justified or not. According 



to calculations in Ref. where the vortices are ignored, 
the presence of an antiferromagnetic order may lead to 
a significant deviation of the d-vector from gk and result 
in some reduction of the paramagnetic depairing effect. 

At present, the best candidate for applying the 
present results in the cubic case will be the familily of 
Li2(Pd3-2;Ptj;)B d, HOl- According to a recent study of 
_ffc2-curves of these materials [2C| . however, they seem 
to be well explained in the weak coupling approximation 
with no paramagnetic effect and in clean limit. Since 
Li2Pd3B is believed to be in the purely s-wave pairing, 
this fact may suggest an extremely weak paramagnetic ef- 
fect in these materials. Nevertheless, the transverse com- 
ponent of the magnetization, stressed in sec. Ill, might be 
measurable, and its experimental search is hoped. 
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